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(54) Fractional sample timing error estimation for W-CDMA 



(57) There is disclosed a method of estimating the 
timing error of a received signal in a spread spectrum 
communication system comprising: cross-correlating a 
portion of the received signal with a known sequence 
thereby generating a sampled cross-correlation func- 
tion; determining the sample having the highest value 
of the sampled cross-correlation function; and estimat- 
ing the timing position of the actual peak of the cross- 
correlation function based on the sample point peak val- 
ue of the sampled cross-correlation function and the val- 
ues of the cross-correlation function at sample points 
before and after the sample generating the peak value. 
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Description 

Field of the Invention 

[0001] The present invention relates to a technique for 
estimating the fractional sample timing error in the re- 
ceiver of a communication system utilising spread-spec- 
trum techniques, and particularly but not exclusively to 
such a technique utilised in a wide-band code division 
multiple access system (W-CDMA). 

Background to the Invention 

[0002] Sampled data W-CDMA receivers sample the 
incoming signal stream from the antenna, after pre- 
processing such as amplification and down conversion, 
at an over-sampling rate equal to a multiple of the nom- 
inal chip rate. If this over-sampling rate is L, then L sam- 
ples per chip are made. In order to de-spread the signal 
every Lth sample is taken and correlated with a locally 
generated spreading sequence that corresponds to the 
known spreading sequence associated with the symbol 
to be received. In conventional receivers, a tracking 
feedback loop is employed that continuously monitors 
the timing error between the locally generated spread- 
ing sequence and the spreading sequence of the re- 
ceived signals. The feedback loop adjusts the local 
spreading sequence generator fine timing to maximise 
the correlation value and hence maximise the signal to 
noise ratio of the de-spread data. Here fine timing refers 
to the fractional chip period timing offset between the 
received signal spreading code and the locally generat- 
ed spreading code replica. Early-late gate code tracking 
loops typify this approach. Examples of this are dis- 
cussed in Chapter 5 of the book "Introduction to Spread 
Spectrum.Communications", Ziemer & Peterson, Pren- 
tice-Hall, 1995. 

[0003] If the over-sampling rate is high enough, i.e. if 
the value of L is high enough, then the timing correction 
can be made by advancing or retarding the local code 
generator by one sample. This technique ensures that 
the timing accuracy is within K» of a chip sample, which 
is equivalent to 1/(2L) of a chip. 
[0004] On the other hand if the value of L is low, say 
4, then the above-described feedback technique will in 
general provide relatively coarse timing correction. This 
is because the correlation peak, which represents the 
optimum sampling point, may be between successive 
samples. The coarseness will introduce losses in corre- 
lation leading to a reduced de-spread data signal-to- 
noise ratio. 

[0005] When low over-sampling rates are employed, 
it is necessary to use some form of interpolation to pro- 
vide an estimate of the signal at the optimum sampling 
point, which lies between samples. Therefore, an esti- 
mate of the fractional sampling delay is required to drive 
an interpolator. 

[0006] A quadratic approximation/interpolation tech- 
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nique has been proposed in "Block Processing Feedfor- 
ward Symbol Tining Estimator for Digital Modems". L.P. 
Sabel and W.G. Cowley. Electronic Letters. 4 th August 
1994, Vol. 30 No. 1 6. pp 1 273-1274. This technique em- 
5 ploys a technique for estimating the peak of a correlation 
function in narrow-band (non-spread) signals. 
[0007] It is an object of the present invention to pro- 
vide a technique for more accurately estimating the tim- 
ing offset of a received spread-spectrum signal, which 

10 technique can provide improved results when the over- 
sampling in the receiver is low. 

Summary of the Invention 

is [0008] According to the invention there is provided a 
method of estimating the timing error of a received sig- 
nal in a spread spectrum communication system com- 
prising: cross-correlating a portion of the received signal 
with a known sequence thereby generating a sampled 

20 cross-correlation function; determining the sample hav- 
ing the highest value of the sampled cross-correlation 
function; and estimating the timing position of the actual 
peak of the cross-correlation function based on the sam- 
ple point peak value of the sampled cross-correlation 

25 function and the values of the cross-correlation function 
at sample points before and after the sample generating 
the peak value. 

[0009] The timing position of the actual peak may be 
approximated by fitting a curve to said sample points. 
30 The sample points before and after the sample gener- 
ating the peak value may be those immediately adjacent 
the peak value. The curve may be a quadratic curve. 
[0010] The step of estimating the timing position of the 
actual peak preferably involves applying the formula: 

35 

3R, -4R 2 + R 3 

40 where t p = the estimated timing position of the peak, f? 7 
= the correlation value at one sample less than the 
peak , R 2 = the correlation value adjacent to the peak, 
and R 3 = the correlation value at one sample more than 
the peak. , 

45 , [0011] The method may further comprise the step of 
adjusting the timing position of the received signal in de- 
pendence on said estimated timing position. The timing 
errors of a plurality of multi-paths of a received signal 
may be estimated, and a plurality of timing position ad- 

so justed received signals may be provided to a respective 
plurality of rake fingers, each rake finger being associ- 
ated with one of the respective plurality of multi-paths. 
[0012] The cross-correlation function may be aver- 
aged for a sequence of received signal portions. 

55 [001 3] In a non-coherent system, the method may fur- 
ther comprise the step of applying a non-linearity func- 
tion to the cross-correlation values. In a coherent sys- 
tem, the method may further comprise the step of mul- 
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tiplying the cross-correlation values with a known com- 
plex sequence. 

[0014] The cross-correlating step may comprise gen- 
erating a continuous cross-correlation function, and 
sampling the continuous cross-correlation function to 
generate the sampled cross-correlation function. 
[0015] According to another aspect of the present in- 
vention there is provided a receiver of a spread spec- 
trum communication system, circuitry for estimating the 
timing error of a received signal, including: cross-corre- 
lation circuitry connected to input the received signal 
and a stored known spreading sequence, and for gen- 
erating a sampled cross-correlation function; peak de- 
tection circuitry, connected to input the sampled cross- 
correlation function, for determining the sample having 
the highest correlation value; and timing estimation cir- 
cuitry, connected to the cross-correlation circuitry and 
the peak detection circuitry, for estimating the timing po- 
sition of the actual peak of the cross-correlation function 
based on the sample point peak value of the sampled 
cross-correlation function and the values of the cross- 
correlation function at sample points before and after the 
sample generating the peak value. 
[0016] The timing estimation circuitry may estimate 
the position of the actual peak by fitting a curve to said 
sample points. The sample points before and after the 
sample generating the peak value are those immediate- 
ly adjacent the peak value. The curve is a quadratic 
curve. 

[0017] The timing estimation circuitry may apply the 
formula: 

" 3R 1 - 4R 2 + R 3 

where tp = the estimated timing position of the peak rel- 
ative to sample location R 1t R 1 = the correlation value 
at one sample less than the peak, R 2 = the correlation 
value adjacent to the peak and which has the maximum 
value, that is R^R'p and R 2 >Ri, and R 3 = the correla- 
tion value at one sample more than the peak. 
[0018] The circuitry may further comprise timing ad- 
justment circuitry connected to receive the received sig- 
nal and to adjust the timing position thereof responsive 
tp the estimated timing position of the actual peak. 
[0019] The benefit of the invention is two-fold. Firstly, 
the fractional timing corrections are made in a feed-for- 
ward manner. In feedback systems, variation of re- 
ceived signal strength corresponds to a variation in the 
feedback loop parameters which may lead to instability 
problems. The feed-forward nature of the apparatus 
eliminates this concern. 

[0020] Secondly, the apparatus may be used to deter- 
mine the timing offsets of several multi-paths from a sin- 
gle set of cross-correlation estimates. This eliminates 
the need for the tracking loops for each multi-path de- 
spreader (known as rake fingers in rake demodulator 
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structures) and consequently leads to reduced imple- 
mentation complexity. 

Brief Description of the Drawings 

5 

[0021] 

Figure 1 illustrates a rake receiver structure accord- 
ing to the present invention for de-spreading multi- 
10 paths in an exemplary W-CDMA system; 

Figure 2 illustrates the main elements of the timing 
error estimation circuitry according to a first pre- 
ferred embodiment of the present invention; 
Figures 3(a) to 3(c) illustrate various stages of the 
15 signal processing performed by the fractional timing 
algorithm implemented in the circuitry of Figure 2; 
Figure 4 illustrates an example squared correlation 
function of a received signal; and 
Figure 5 illustrates the main elements of the timing 
error estimation circuitry according to a second pre- 
ferred embodiment of the present invention. 

Description of Preferred Embodiment 

[0022] The invention is described herein with refer- 
ence to an exemplary embodiment. One skilled in the 
art will appreciate, however, that the invention is not lim- 
ited to such an embodiment, and variations may apply 
in different implementations. 

[0023] Figure 1 illustrates a portion of a receiver in a 
spread-spectrum W-CDMA system. The signal oh line 
2 is provided to the receiver circuitry by a front-end proc- 
essor of the receiver, which down-converts the received 
signal from RF to nominal baseband, then converts the 
nominal baseband analogue signal into digital format by 
sampling, and chip-matched filters the received signal. 
The signal on line 2 is thus a sampled, chip-matched 
filtered version of the wanted signal received at the front- 
end of the receiver. The rate at which the received signal 
is sampled is nominally a multiple of the chip rate of the 
received spread-spectrum signal. 
[0024] The portion of the receiver circuitry shown in 
Figure 1 is the circuitry for recovering transmitted sym- 
bols. Generally designated by reference numeral 10 is 
a conventional rake receiver comprising a plurality 8a to 
8 N of multi-path de-spreading circuits ("rake fingers"). 
The transmitted signal for a particular caller travels from 
the transmitter to the receiver in the air interface via a 
possible plurality of paths. These paths may include a 
direct path from transmitter to receiver, and also include 
signals reflected along the path between transmitter and 
receiver. Each receiver is typically provided with an ap- 
propriate number of rake receivers 8a to 8 N to be able 
to decode a sufficient number, N, of multi-paths for a 
particular caller, to recover a strong replication of the 
transmitted signal at the receiver. 
[0025] According to the invention, the conventional 
rake receiver is further adapted and provided with a tim- 
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ing error estimation circuitry 4 and a plurality of multi- 
path interpolators 6a to 6 N . each one of the multi-path 
interpolators being associated with a respective one of 
the rake fingers 8a to 8 N . 

[0026] The sampled received signal is provided as an 
input on line 2 to each of the multi-path interpolators 6a 
to 6 N , and to the timing error estimation circuitry 4. The 
timing error estimation circuitry 4 provides an N plurality 
of outputs on lines 12, one to each of the multi-path in- 
terpolators 6a to 6 N , on a respective signal line 12a to 
1 2 N . The timing error estimation circuitry 4 also provides 
an N plurality of outputs on lines 13. one to each of the 
multi-path de-spreading circuits 8a to 8 N , on a respec- 
tive signal line 13a to 13 N . 

[0027] As will be described further hereinbelow, the 
timing error estimation circuitry 4 estimates both frac- 
tional and integer timing errors of multi-paths. The plu- 
rality of outputs on lines 1 2 from the timing error estima- 
tion circuitry each indicate a fractional timing error of a 
multi-path (a fractional timing error being less than a 
chip sample period). The plurality of outputs on lines 13 
from the timing error estimation circuitry each indicate 
an integer chip sample period timing error. 
[0028] The multi-path dc-sprcading circuits 8a to 8 N 
each receive as an input the outputs of a respective one 
of the multi-path interpolators 6a to 6 N on lines 16a to 
16 N respectively. The multi-path de-spreading circuitry 
8a to 8 N each present on their respective outputs 1 4a 
to 14 N a de-spread and time-shifted recovered symbol 
estimate of the transmitted symbol. The recovered sym- 
bols on lines 14a to 14 N are then combined in the com- 
biner 18 which presents the recovered final symbol on 
line 20 for further processing in the receiver. 
[0029] *• The main functional elements of the timing er- 
ror estimation circuitry according to the invention are 
shown in Figure 2, for an example of a non-coherent 
system. The timing error estimation circuitry comprises 
cross-correlation circuitry 22, non-linearity circuitry 24, 
averaging circuitry 26, peak location circuitry 28, and al- 
gorithm circuitry 30. 

[0030] The operation of the circuitry to estimate the 
error in the timing of the received signal will now be de- 
scribed with reference to Figures 1 and 2 and also with 
reference to Figure 3, which illustrates the resulting sig- 
nals at various points in the processing stages of the 
timing error estimation circuitry 4. 
[0031] The cross-correlation circuitry 22 receives on 
signal line 2 a sample of the received signal including 
the spreading code incorporated in the signal at the 
transmitter. The cross-correlation circuitry is provided, 
under local control (not shown) with a locally generated 
version of the spreading code used in the transmitter. 
The receiver knows the spreading codes used by the 
transmitter, and provides the appropriate spreading se- 
quence to the cross-correlation circuitry. Only one 
spreading code need be provided to the fractional timing 
circuitry for each user: for a particular caller or user the 
spreading code is the same for all multi-paths, only shift- 



ed in time. 

[0032] In a conventional manner, the cross-correla- 
tion circuitry correlates the spreading code of the re- 
ceived signal with the locally generated spreading code. 
5 The estimated timing offset of the received signal is the 
position where the cross-correlation is a maximum, 
which is generally not coincident with a cross -correlation 
sample value. 

[0033] The thus calculated correlation values are out- 
10 put on lines 32 to the non-linearity circuitry 24. The non- 
linearity circuitry may implement any one of a number 
of functions, including a square law or absolute value 
function. In this example there is described an imple- 
mentation using a square law device or squaring circuit- 
's ry. The invention, however, is not limited to such an im- 
plementation. Thus, the squaring circuitry squares the 
values of the correlation function, and outputs the 
squared values to the averaging circuitry 26 on signal 
lines 34. 

20 [0034] The correlation values are squared in the ex- 
ample of a non-coherent system because the correla- 
tion components are complex. By squaring, the sum of 
the square of the in-phase components and the sum of 
the square of quadrature components is obtained to give 

25 the magnitude of the correlation values. 

[0035] Referring to Figure 4, there is illustrated an ex- 
ample of the squared cross-correlation function which 
is generated by the cross-correlation circuitry 22, and 
non-linearity circuitry 24 on line 34. For each received 

30 sample the cross-correlation circuitry 22 performs a 
number of cross-correlations for different estimates of 
the timing position of the signal. For example, in one 
preferable implementation the cross-correlation circuit- 
ry 22 performs a total of 320 cross -correlations, for 320 

35 estimates of the starting position of the transmitted sym- 
bol in the received sample. Thus, for each received sam- 
ple the cross-correlation circuitry performs a set of 320 
cross-correlations. 

[0036] As can be seen from Figure 4, the generated 
40 squared cross-correlation function has a number of 
peaks, which correspond to the multi-paths associated 
with a particular transmitted symbol. Thus the peak MP1 
corresponds to multi-path 1 , the peak MP2 corresponds 
to multi-path 2, and the peak MP3 corresponds to multi- 
45 path 3. 

[0037] The averaging circuitry 26 averages succes- 
sive sets of correlation values. For example, for each 
successive set of received samples the averaging cir- 
cuitry updates an average correlation function, which 

so gives a representation of the averaged correlation of 
successive received sets of samples over time. 
[0038] The averaging circuitry outputs on signal lines 
36 the average of the squared correlation functions over 
time. It will be understood that the cross-correlation cir- 

55 cuitry 22, the squaring circuitry 24 and the averaging cir- 
cuitry 26 may in fact form a single functional block 21 , 
which outputs the averaged square correlation function; 
where the number of averaged, squared, cross-correla- 



BNSDOCID: <EP 1089452A1 J_> 



7 



EP 1 089 452 A1 



8 



tion functions may be varied, to obtain the required tim- 
ing estimate variance given the received signal-to-noise 
ratio. 

[0039] The thus generated averaged squared cross- 
correlation function has a shape as shown in Figure 3 5 
(a) for one multi-path. It will be appreciated that the av- 
eraged squared cross-correlation function also includes 
the peaks associated with the other multi-paths associ- 
ated with the function of Figure 4, but only one peak is 
shown throughout Figure 3 for ease of explanation. The 10 
time axis in Figure 3(a), as in Figure 4, does not repre- 
sent real time. Rather the time axis represents the time 
relative to an estimated start location, or bit position, of 
the spreading code of the received signal. 
[0040] The thus generated averaged squared cross- »5 
correlation function on signal lines 36 forms an input to 
the peak location circuitry 28. The peak location circuitry 
28 samples the averaged squared cross-correlation 
function at a sample rale having an over-sampling factor 
of L. For coherent detection an over-sampling factor as 20 
low as L=2 may be chosen. For this example of a non- 
coherent system, an over-sampling factor of L=4 is pref- 
erably used. 

[0041 ] That is the sample rate of the cross-correlation 
function is L times the nominal chip rate. The optimum 25 
sampling position is the peak of the cross-correlation. 
The peak of the cross-correlation function is the best 
estimate of where the spreading code is located in the 
received sample set 

[0042] As will be understood by one skilled in the art, 30 
the peak location circuitry 28 will detect the peaks to 
conform with certain criteria. For example, the peak val- 
ue will be determined as being associated with a multi- 
path only if it exceeds a certain threshold value, and the 
peaks of multi-paths have to be a certain distance apart. 35 
The implementation of such peak location circuitry for 
identifying the locations of multi-paths is well-known, 
and is outside the scope of the present invention. In the 
example of Figure 4, the cross-correlation circuitry iden- 
tifies three multi-path locations for a particular transmit- 40 
ted symbol: MPt , MP2, MP3. 

[0043] The averaged squared cross correlation func- 
tion of Figure 3(a) is sampled at times Aj. As can be seen 
in Figure 3(b), this results, when extrapolated, in a func- 
tion where a peak cross-correlation value can be iden- 45 
tif ied for a particular sample. It can be seen that the peak 
value occurs at a sample time a 2 , and this can be de- 
tected by the peak location circuitry in the conventional 
way by detecting a lower cross-correlation value at the 
sample times A-, and a 3 either side of the sample A 2 . 50 
The cross-correlation value at sample time A n is R1, at 
sample time A 2 is R2 , and at sample time A 3 is R3. The 
estimate of the location of the peak for the multi-path 
shown in Figure 3, and the other multi-paths of Figure 
4, is used by the peak location circuitry to generate tim- 55 
ing offset signals associated with each of the multi- 
paths. These signals are generated on signal lines 13a 
to 13 N . Thus the signals on lines 13a to 13 N indicate the 



timing offset to be applied to the respective multi-paths 
so that each of the code generators in the multi-path de- 
spreading circuitry 8a to 8 N can adjust the position at 
which the local spreading code is applied for cross-cor- 
relation. 

[0044] Referring to Figures 3(a) and 3(b) it can be 
clearly seen that the sample time A 2 is in fact to the right 
of where the actual peak of the averaged squared cross- 
correlation function is. Clearly, the optimum sampling 
position lies between the sample times A 1 and A 2 . It is 
assumed that this 'coarse' timing acquisition positions 
the multi-path peaks to within half a chip period of the 
actual peak. 

[0045] To achieve the best fractional timing error cor- 
rection, it is necessary to make a good estimation of the 
time difference between the time of the actual peak of 
the cross-correlation function as designated by time t p 
in Figure 3(c), and the time A 2 . 

[0046] For systems where a high over-sampling rate 
is used the location of the sample A 2 will approximate 
more closely the actual peak. Where a low over-sam- 
pling rate is used, the difference between the timing po- 
sition of the actual peak and the timing of the sample A 2 
may introduce a timing error which significantly de- 
grades system performance. 

[0047] To overcome this error in the estimate of the 
location of the peak of the cross-correlation function in 
this preferred embodiment, the peak location detector 
circuitry outputs the peak values R1 , R2 and R3 on lines 
38 to an algorithm circuitry 30. The object of the algo- 
rithm circuitry 30 is to estimate the position of the peak 
to a fraction of a sample. The algorithm^ rcuitry as- 
sumes that the three values R1 , R2 and R3 are points 
on ah underlying smooth curve, i.e. the algorithm circuit- 
ry works by assuming that the three points, taken at 
times A v a 2 and A 3 , are samples of a smoother under- 
lying curve. 

[0048] Thus the algorithm circuitry 30 fits a polynomial 
to these three points. In a preferred embodiment, a pa- 
rabola can be calculated which fits the points exactly 
without recourse to minimum mean squared error tech- 
niques. The parabola is shown in Figure 3(c), and isgen- 
erally designated by reference numeral 40. 
[0049] The parabola joining the three points R1, R2, 
and R3 can be defined by the general equation at 2 + bt 
+ c = R , where t is the sample time and R is the aver- 
aged squared cross correlation output, i.e. a particular 
value of the averaged squared cross-correlation func- 
tion at the particular sample time. Given the received 
cross-correlation values R h where in this case *={1 .2,3}, 
it is possible to form the following simultaneous equa- 
tions: 

aA* + f?A 1 + c= fl, 

2 

aA 2 + b\ 2 + c — R 2 
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aA 3 + bA 3 + c = R 3 

[0050] Inserting the values 0.1 and 2 for A v A 2 and 
A 3 results in the following matrix equation: 



"0 
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b 
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c 







[0051 ] Inverting the three by three matrix gives the so- 
lution for a, b, and c. The solution is found from: 



''/: 
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2 


-A 
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0 






c* 



[0052] The peak of the parabola is at the turning point, 
found by consideration of the first derivative of the par- 
abolic equation, and given by: 

t p = -b/2a 

[0053] Substituting for b and a gives: 
3R, - 4R 2 + R 3 

{ p = a, -2P 2 + r 3 

[0054] The algorithm circuitry then calculates the tim- 
ing difference between the determined value t p and the 
first sample value A, to give the calculated fractional tim- 
ing error. 

[0055] While other similar formulae may be derived 
by using sample locations R 2 and R 3 as the reference 
point, the above formula is used for example purposes. 
[0056] In this example the integer offset used on lines 
13 corresponds to the location of the maximum value of 
the averaged squared cross-correlation, i.e. A 2 , where- 
as the fractional part is relative to A 1 . Hence in this case 
the estimate tp must be adjusted by the subtraction of 
a sampling period. Alternatively, the above equation 
may be formulated using A 2 , as the reference point rath- 
er than A -) as used in the above derivation. 
[0057] This timing error is calculated for each of the 
peaks of the averaged squared cross-correlation func- 
tion, only one of which is shown in Figure 3(a), thus gen- 
erating N fractional error timing estimates associated 
with each of the N multi-paths. This timing error infor- 
mation is then generated on each of the signal lines 1 2a 
to 12N and sent to the respective multi-path interpola- 
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tors 6a to 6 N . 

[0058] It should be noted that the peak location t p is 
still an estimate: it is based on the assumption that each 
of the peaks of the average squared correlation function 

5 has a parabola shape. However, it gives a much more 
accurate estimate than the peak sample value. 
[0059] The respective interpolators 6a to 6 N use the 
timing error information provided on lines 12a to 12 N to 
shift the signal on line 2 by the appropriate amount (ei- 

10 ther forward or backward) prior to input to one of the 
fingers of the rake receiver. 

[0060] Thus the timing error estimation circuitry 4 in 
combination with the multi-path interpolators 6a to 6 N 
provide fractional sample timing corrections. The timing 
15 corrections are employed in a feed-forwards basis as 
shown in Figure 1 . In this structure the fractional timings 
are passed to respective interpolators which produce an 
output corrected by the fractional sample error. 
[0061] In the above description the inventive tech- 
no nique has been described with reference to an example 
of non-coherent system. However, the invention is more 
generally applicable, and one skilled in the art will ap- 
preciate that the invention is not limited to such a sys- 
tem. In an alternative implementation described herein - 
25 after, the technique of the present invention is described 
applied by way of further example in a preferred imple- 
mentation of a coherent system. 
[0062] In a coherent system, known data (e.g. pream- 
ble or pilot data symbols) is used which provides im- 
30 proved performance (lower timing offset estimate vari- 
ance). In the coherent system; instead of calculating the 
squared magnitudes of the cross-correlation values of 
the received signal as described in the above example, 
the complex product of the cross-correlation values and 
35 the complex conjugate of the known data symbols of the 
received signal is calculated. 

[0063] For completeness, referring to Figure 5, the 
main elements of the implementation of the timing error 
estimation circuitry 4 for a coherent system are shown. 

40 Like reference numerals are used to indicate elements 
which are the same as those shown in Figure 2. The 
only difference is in the functional block 21, where an 
additional block for multiplying by a known symbol se- 
quence 50 is included, in addition the averaging circuitry 

45 and non-linearity circuitry 24 swap positions. The multi- 
plying operation in block 50 is a complex multiplication 
step. 

[0064] In the above-described example, two particu- 
lar implementations of the invention for coherent and 
so non-coherent detection are described. One skilled in the 
art will appreciate that variations may be made to the 
described implementations without departing from the 
scope of the invention. 

[0065] In the preferred implementation, where the 
55 peak estimation is based on the fitting of the detected 
points to the equation of the above-stated parabola, a 
sample rate of four samples per chip is required at the 
input of the peak location circuitry 28. However, in a pre- 
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ferred implementation, the received signal on line 2 has 
a sample rate of only two samples per chip. In such a 
scenario, prior to the cross-correlation circuitry 22 an in- 
terpolation stage may be introduced such that the signal 
on line 2 is interpolated by a factor of 2 to provide a signal 5 
to the cross-correlation circuitry 22 having a sample rate 
of four samples per chip. 

[0066] Alternatively, such interpolator having an inter- 
polation factor of 2 may be introduced after the cross- 
correlation circuitry 22. This is advantageous, since it w 
enables the cross-correlation circuitry 22, which is pref- 
erably implemented using a code-matched filter, to be 
half, the size it would be if the interpolation was per- 
formed at its input. 

[0067] in a coherent system, the interpolation by a '5 
factor of two could be implemented at any point prior to 
the peak location circuitry 28, since the process is linear 
throughout. In a non-coherent system, the interpolation 
by a factor of two must be implemented prior to the non- 
linearity circuitry 24. 20 
[0068] In addition, although hereinabove a particular 
set of equations is described for approximating a curve 
to the three detected sample points, other variations 
may exist. Any technique for approximating a curve to 
a set of points may be used, and then formula for de- 25 
tecting the peak of that curve applied. 
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1 . A method of estimating the timing error of a received 
signal in a spread spectrum communication system 
comprising: cross-correlating a portion of the re- 
ceived signal with a known sequence thereby gen- 
erating a sampled cross-correlation function; 35 
determining the sample having the highest value of 
the sampled cross-correlation function; and 
estimating the timing position of the actual peak of 
the cross-correlation function based on the sample 
point peak value of the sampled cross-correlation 40 
function and the values of the cross-correlation 
function at sample points before and after the sam- 
ple generating the peak value. 

2. The method of claim 1 wherein the timing position 45 
of the actual peak is approximated by fitting a curve 

to said sample points. 

3. The method of claim 1 or claim 2 wherein the sam- 
ple points before and after the sample generating so 
the peak value are those immediately adjacent the 
peak value. 

4. The method of claim 2 or claim 3 wherein said curve 

is a quadratic curve. 55 

5. The method of any one of claims 1 to 4, wherein the 
step of estimating the timing position of the actual 



peak involves applying the formula: 



where t p - the estimated timing position of the peak, 
R 1 = the correlation value at one sample less than 
the peak , R 2 = the correlation value at the peak, 
and R 3 =[he correlation value at one sample more 
than the peak. 

The method of any one of claims 1 to 5 further com- 
prising the step of adjusting the timing position of 
the received signal in dependence on said estimat- 
ed timing position. 

The method of claim 6 wherein the timing errors of 
a plurality of multi-paths of a received signal are es- 
timated, and a plurality of timing position adjusted 
received signals are provided to a respective plu- 
rality of rake fingers, each rake finger being asso- 
ciated with one of the respective plurality of multi- 
paths. 

The method of any one of claims 1 to 7 in which the 
cross-correlation function is averaged for a se- 
quence of received signal portions. 

The method of any one of claims 1 to 8 in a non- 
coherent system, further comprising the step of ap- 
plying a non-linearity function to the cross-correla- 
tion values. 



10. The method of any one of claims 1 to 8 in a coherent 
system, further comprising the step of multiplying 
the cross-correlation values with a known complex 
sequence, 

11. The method of any one of claims 1 to 10, wherein 
the cross-correlating step comprises generating a 
continuous cross-correlation function, and sam- 
pling the continuous cross-correlation function to 
generate the sampled cross-correlation function. 

12. A method of recovering a transmitted symbol in a 
received signal in the receiver of a spread-spectrum 
system including the steps of estimating the timing 
error of the received signal according to any one of 
claims 1 to 11, adjusting the timing position of the 
received signal in accordance with the determined 
error and de-spreading the time-adjusted received 
signal. 

13. In a receiver of a spread spectrum communication 
system, circuitry for estimating the timing error of a 
received signal, including: cross-correlation circuit- 
ry connected to input the received signal and a 
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stored known spreading sequence, and for gener- 
ating a sampled cross-correlation function peak de- 
tection circuitry, connected, to input the sampled 
cross-correlation function, for determining the sam- 
ple having the highest correlation value; and timing 
estimation circuitry, connected to the cross-correla- 
tion circuitry and the peak detection circuitry, for es- 
timating the timing position of the actual peak of the 
cross-correlation function based on the sample 
point peak value of the sampled cross-correlation 
function and the values of the cross-correlation 
function at sample points before and after the sam- 
ple generating the peak value. 

14. The circuitry of claim 13 wherein the timing estima- 
tion circuitry estimates the position of the actual 
peak by fitting a curve to said sample points. 

15. The circuitry of claim 13 or claim 14 wherein the 
sample points before and after the sample generat- 
ing the peak value are those immediately adjacent 
the peak value. 

16. The circuitry of claim 14 or 1 5 wherein said curve is 
a quadratic curve. 

17. The circuitry of any one of claims 13 to 16, wherein 
the timing estimation circuitry applies the formula: 



3R r 4R 2 +F? 3 



where t p = the estimated timing position of the peak, 
R 1 = the correlation value at one sample less than 35 
the peak , R 2 = the correlation value at the peak, 
and R 3 = the correlation value at one sample more 
than the peak. 
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18. The circuitry of any one of claims 13 to 17 further 
comprising timing adjustment circuitry connected to 
receive the received signal and to adjust the timing 
position thereof responsive to the estimated timing 
position of the actual peak. 
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